Reconnection rate for the steady-state Petschek model 
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Abstract 

Reconnection rate is found for the canonical simplest case of steady- 
state two-dimensional symmetric reconnection in an incompressible plasma 
by matching of outer Petschek solution and internal diffusion region solution. 
The reconnection rate obtained naturally incorporates both Sweet-Parker 
and Petschek regimes, the latter seems to be possible only for the case with 
strongly localized resistivity. 



I. INTRODUCTION 

Magnetic reconnection is an energy conversion process which occurs in astro- 
physical, solar, space and laboratory plasmas (e.g., Hones 1 ; Priest 2 ). First attempts 
to explain the fast energy release in solar flares based on pure resistive magnetic field 
dissipation (Sweet 3 ; Parker 4 ) showed that the energy conversion rate is estimated 
as 1/ V-Rem, where 



is the global Reynolds number, L is the half-length of reconnection layer, Va is 
Alfvenic velocity, and 77 is resistivity. For typical conditions in the solar corona the 
Sweet-Parker rate turns out to be orders of magnitudes too small when compared 
to experimental data. 

In 1964 Petschek 5 pointed out that in a highly-conducting plasma dissipation 
needs only to be present within a small region known as the diffusion region, and 
energy conversion occurs primarily across non-linear waves, or shocks. This gives 
another estimation of the maximum reconnection rate l/\nRe m which is much more 
favorable for energy conversion. 

Unfortunately, up to the present it is still unclear which conditions make 
Petschek-type reconnection to be possible and which are responsible for the Sweet- 
Parker regime. The fact is that numerical simulations (Biskamp,1986, Scholer, 1989) 
were not able to reproduce solution of Petschek type but rather were in favor of 
Sweet-Parker solution unless the resistivity was localized in a small region (e.g., Sc- 
holer 1989, Yan, Lee and Priest, 1992, Ugai, 1999). The laboratory experiments also 
seem to observe Sweet-Parker regime of reconnection (Uzdensky et al., 1996, Ji et 
al.,1999). 

From the mathematical point of view the problem of reconnection rate is con- 
nected with the matching of a solution for the diffusion region where dissipation is 
important, and solution for the convective zone where ideal MHD equations can be 
used. But up to now this question is still not resolved even for the canonical simplest 
case of steady-state two-dimensional symmetric reconnection in an incompressible 
plasma. 

It is the aim of this paper to present a matching procedure for the canonical 
reconnection problem. The reconnection rate obtained from the matching turns out 
to incorporate naturally both Petschek and Sweet-Parker regimes as limiting cases. 



Petschek solution 

We consider the simplest theoretical system consisting of a two-dimensional 
current sheet which separates two uniform and identical plasmas with oppositely 
oriented magnetic fields ±B . Petschek (1964) pointed out that the diffusion region 
can be considerably smaller than the whole size of the reconnection layer and that 
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the outer region contain two pairs of standing slow shocks. These shocks deflect 
and accelerate the incoming plasma from the inflow region into two exit jets wedged 
between the shocks (see Figure 1). This jet area between the shocks with accelerated 
plasma is traditionally called outflow region. 

In the dimensionless form the Petschek solution can be presented as follows 
(Petschek, 1964, for details see Vasyliunas, 1975): 

Inflow region: 

v x = 0; v y = -e, (2) 

As 1 As x 

B x = l In = ; B y = — arctan-. (3) 

7i \fx 2 + y 2 7i y 

Outflow region: 

v x = 1; v y = 0; B x = 0; B y = e. (4) 
Equation of shock in the first quadrant is the following: 

y = ex. (5) 

Here x, y are directed along the current sheet and in the perpendicular direction, 
respectively. We normalized the magnetic field to B , length to L, plasma velocity 
to Alfvenic velocity Va, and electric field E to Alfvenic electric field Ea = VaBq. 

The reconnection rate 

e = E/E A « 1 (6) 

is supposed to be a small parameter of the problem. 

Expressions (0-^) are the asymptotic solution with respect to e of the MHD 
system of equations 

(v • V)v = -VP + (B • V)B, (7) 

E+tvxBH^curlB, (8) 

divB = 0, divv = 0, (9) 

and the Rankine-Hugoniot shock relations in the limit Re m — > oo. Petschek did not 
obtain a solution in the diffusion region, instead he estimated maximum reconnection 
rate as l/\nRe m of using some simple physical suggestion . Generally speaking, this 
implies that the Petschek model gives any reconnection rate from Sweet-Parker value 
1 / V Re m up to l/lni?e m , and it is still unclear whether Petschek reconnection faster 
than Sweet-Parker reconnection is possible. The problem can be solved by matching 
of a solution for the diffusion region and Petschek solution 



Diffusion region scaling 

We renormalize the MHD equations to the new scales B' , V A , E' A = B' V A , where 
all quantities are supposed to be taken at the diffusion region upper boundary, and 
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at the half length of the diffusion region 1$. We have to use the dissipative MHD 
equations (0^0) for the diffusion region with Reynolds number 

Re' m = (10) 
V 

and electric field E = e' . 

The scaling for the diffusion region is similar to that for the Prandtl viscous layer 
(see Landau and Lifshitz, 1985): 

x', B' x , <, P'~0(1), 



y', B' y , v' y , e' ~ l/y/Re' m . (11) 
Consequently, the new boundary layer variables are the following: 



y = y'^Re'm, B y = B'^Re' m , v y = v' y ^Re' rn , e = e'^Re' m . (12) 
The diffusion region Reynolds number is supposed to be Re' m >> 1, and therefore in 



the zero-order with respect to the parameter 1/ J Re' m the boundary layer equations 
turn out to be: 

y dv^_ I ~ dv^_ _ g 8B X _ g dB x _ dP(x) 
x dx V dy x dx V dy dx ' \ / 

divB = 0, divv = 0, (14) 
P = P(x), (15) 
v y B x -v x B y -e = fi{x,y)^, (16) 

where fj(x,y) is the normalized resistivity of the plasma with maximum value 1. 

Unfortunately, the appropriate exact solutions of the boundary layer equations 
(|13 HT6|) are unknown, therefore we have to solve the problem numerically. The main 
difficulty is that the internal reconnection rate e is unknown in advance and has 
to be determined for given resistivity fj(x,y), given total pressure P(x), and B x (x) 
given at the upper boundary of the diffusion region. In addition, the solution must 
have Petschek-type asymptotic behaviour (0-§[) outside of the diffusion region. 

Although we are looking for a steady-state solution, from the simulation point of 
view it is advantageous to use relaxation method and solve numerically the following 
unstationary system of boundary layer MHD equations: 



§v _l_ i) §0*. _l_ 7, Hik _ R 9B * — R — _ dP ( x ) ( 1 7^ 

dt ~r Ux dx ^ U V dy JDx dx D V dy ~ dx ' \ Ll ) 

^ = curl(v x B) — curl (f](x, y) curlB^ , (18) 

divB = 0, divv = 0. (19) 
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Starting with an initial MHD configuration under fixed boundary conditions we look 
for convergence of the time-dependent solutions to a steady state. 

As initial configuration we choose a X-type flow and magnetic field: v x = x , v y = 
—y, B x = y, B y = —x. The distribution of the resistivity is traditional (see 
(Ugai,1999, Scholer 1985)): 

rj( x , y) = de ( - s ^ 2 - s ^ 2) + f, (20) 

with d + / = 1 where coefficient d describes inhomogeneous resistivity, and / is 
responsible for the background resistivity. 

The problem under consideration consists essentially of two coupled physical 
processes: diffusion and wave propagation. To model these processes, two-step with 
respect to time numerical scheme has been used. At first, convectional terms were 
calculated using the Godunov characteristic method, and then the elliptical part 
was treated implicitly. 

Calculations were carried out on a rectangular uniform grid 100 x 145 in the first 
quadrant with the following boundary conditions: 

Lower boundary: symmetry conditions dv x /dy = 0, v y — 0, B x = 0; induction 

equation ( P^|) has been used to compute the B y component at the x-axis. 

Left boundary: symmetry conditions v x = 0, dv y /dx = 0, dB x /dx = 0, B y = 0. 

Right boundary: free conditions dv x /dx = 0, dv y /dx = 0. 

Upper (inflow) boundary: v x = 0, B x — 1. 

Note, that this implies that we do not prescribe the incoming velocity, and hence 
the reconnection rate: the system itself has to determine how fast it wants to recon- 
nect. 

The total pressure can be fixed to 1 in the zero-order approximation: P — 1. 

Let us discuss the result of our simulations. For the case of localized resistivity 
where we chose d = 0.95, / = 0.05, s x = s y = 1 in the equation fl2"Dp, the 
system reaches Petschek steady state (see Figure 2) with clear asymptotic behaviour, 
pronounced slow shock, and the reconnection rate turns out to be e ~ 0.7. 

From the other hand, for the case of homogeneous resistivity d = 0, / = 1, the 
system reaches Sweet-Parker state (see Figure 3) with much less reconnection rate 
i ~ .25 even if the Petschek solution has been used as initial configuration (see also 
(Ugai, 1999, Scholer, 1989)). This seems to imply that Petschek-type reconnection is 
possible only if the resistivity of the plasma is localized in a small region, and for 
constant resistivity the Sweet-Parker regime is realized. 

The size of the diffusion region can be defined as the size of the region where 
the convective electric field E = v x B (which is zero at the origin) reaches the 
asymptotic value e (or, some level, say 0.95e). For the case of localized resistivity 
Id practically coincides with the scale of the inhomogeneity of the conductivity. In 
principal, there might be a possibility to produce Petschek-type reconnection with 
constant resistivity using a highly inhomogeneous behaviour of the MHD parameters 
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at the upper boundary (narrow stream, for example, see Chen et al.,1999), and then 
Id has the meaning of the scale of this shearing flow or other boundary factor which 
causes the reconnection. 



Matching procedure 

We have only a numerical solution for the diffusion region, and this makes it dif- 
ficult for the matching procedure because the latter needs an analytical presentation 
of the solutions to be matched. The only way out left is to continue the diffusion 
region solution to the inflow region using dates known from the simulation distribu- 
tion of the By component along the upper boundary of the diffusion region. Then 
try to match the solutions in the current free inflow region at the distance r ~ Id 
(see Figure 1). 

As can be seen from equation ([!]) the B x component of the Petschek solution 
diverges at the origin B x — > — oo when r = \/x 2 + y 2 — > 0. This singularity is a 
consequence of the fact that dissipation actually has not been taken into account for 
the solution (0- 15[) which is nevertheless still valid untill the distances of the order 
of the size of diffusion region is Id- 

In order to be adjusted to the Petschek solution, the B' y component must have 
the following limit for x/ld — > oo at the upper boundary of the diffusion region : 

B' y {x/l d ) -> 2e. (21) 

We can obtain the asymptotic behaviour of B' x for r > Id region using a Poisson-like 
integral presentation: 



,,, n , 1 /■+*> dB'yHx,0) 
B> x (x f ,y f ) = B>+- V M n 

71 J-oo OX i 



oodB'W(x,0) J(x-x) 2 + y 2 

1 dx 



d 



+ i /+ ~ 8Biy, o) r y^Tg + g-jgi = 

7r J-oo at, I l d x 2 + y 2 J 

# + -lnf + 0(l/r), (22) 

7T Id 

where ^ = x/ld- This gives an outer expansion for the inner solution. On the other 
hand a convective solution (§) can be rewritten in the following form in order to 
determine the inner expansion of the outer solution: 

4e, I As L Ae l d , s 

B x = l In — = 1 In- In-. 23 

7r r n Id 71 r 

Equating these two asymptotic expansions we obtain the matching relation: 

Ae L 

5=1- -In-, (24) 

7T Id 
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Now everything is ready to determine the reconnection rate. The electric field must 
be constant in the whole inflow region, hence 



v'B' Q = vB , (25) 
e'B'J = sBl (26) 

where the definition of the reconnection rates e' = v'/B' Q , e = v/Bq has been used. 
Bearing in mind that e' = ej VR~e' m (see scaling (|l|) we obtain: 



s~Bf 2 = eB^Y^. (27) 

Substituting B' from the equation (^4j) we determine finally the following equa- 
tion for the reconnection rate e : 



s(l--\n^f 2 = sJReJf, (28) 
n l d \ L 

where Re m is the global Reynolds number (TO), and the internal reconnection rate e 
has to be found from the simulation of the diffusion region problem. 
For small e there is an analytical expression: 

(29) 



ReJ-i + ^elnf 



Here e is an internal reconnection rate determined from the numerical solution: 
e ~ 0.7. 



Discussion and conclusion 

Equations (^8],^9]) give the unique reconnection rate for known parameters of 
the current sheet L, B , Va, i], Id- For sufficiently long diffusion region such 
that \J Re m l -f^ » ^edn^, the equation fl29"|) corresponds to Sweet-Parker regime 

e ~ ej J Re m l -f;. For the opposite case of resistivity constrained in a small region e ~ 
|/ In j- we have Petschek reconnection. Hence, reconnection rate (|28| , |29| ) naturally 
incorporates both regimes obtained in simulations (Scholer, Ugai, Biskump). 

We were not able to reproduce Petschek regime using variation of MHD parame- 
ters at the upper boundary with homogeneous resistivity, a probably solution (Chen, 
1999) of this problem either is essentially time-dependent or corresponds to the case 
of strong reconnection. According to our simulations, for Petschek state to exist a 
strongly localized resistivity is needed, and for the spatially homogeneous resistivity 
Id = L Sweet-Parker regime seems to be always the case. This result resolves old 
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question about conditions which are necessary for Petschek-type reconnection to 
appear. 

It is interesting that for the deriving of equations (|28| , |29D the only value which has 
been actually used is the internal reconnection rate e obtained from the numerical 
solution, but the distribution of the B y component along the upper boundary of the 
diffusion region does not contribute at all (besides asymptotic behaviour (|22|)) in the 
zero-order approximation considered above. Of course, from the mathematical point 
of view it is important that diffusion region solution exists and has Petschek-like 
asymptotic behaviour (fJ-§). 

The strongly localized resistivity is often the relevant case in space plasma ap- 
plications, but for the laboratory experiments where the size of a device is relatively 
small the Sweet -Parker regime is expected. 
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Figure Captions 

Figure 1: Scheme of matching of the outer Petschek solution and diffusion region 
solution. 

Figure 2: Configuration of magnetic field lines (solid line) and stream lines 

(dashed line) for the numerical simulation of the diffusion region. 
Figure 3: Three-dimensional plot of current density shows Petschek shock 
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